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Abstract 

We give a method for computing the ensemble average of multiplicative class func- 
tions over the Gaussian ensemble of real asymmetric matrices. These averages are 
expressed in terms of the PfafBan of Gram-like antisymmetric matrices formed with 
respect to a skew-symmetric inner product related to the class function. 

1 Introduction 

In the 1960's Ginibrc introduced three statistical ensembles of matrices whose entries are 
chosen independently with Gaussian probability from (resp.) M, C, and Hamilton's Quater- 
nions [4]. These ensembles are respectively labeled GinOE, GinlJE and GinSE in analogy 
with their Hermitian counterparts. Ginibre introduced a physical analogy between GinUE 
and two-dimensional electrostatics, but gave no applications for the other two ensembles. 
Since their introduction, many applications have been found for GinOE and GinSE. 

Here we report a method for determining ensemble averages over GinOE of certain 
functions which are constant on similarity (conjugacy) classes. GinOE is the space oi N x N 
real matrices R^^^ together with the probability measure v given by 

Js 

where fi is Lebesgue measure on R^^^, and (27r)^'/2. Our goal is to find 

(*>:=23jv' J *(X)cxp(-Tr(XTx)/2)d/x(X), 

RNXN 

where * : R^^^ ^ R is (i) constant on similarity classes: ^i{AXA^''-) = for all 

invertible A G R^^^, and (ii) there exists a function : C ^ R such that if D is a 
diagonal matrix with entries 71, 72, ... , then '^{D) — ip ill) 4' (12) ■ ■ • V'(7Af)- We remark 
that /z-almost every matrix in R^^^ is similar to a diagonal matrix and consequently 
is uniquely determined by ip. 

As an example, when ^^(7) = 7": we are in the well-studied situation 'i'(X) — (detX)" 
[8][1]. 

Ginibre's original interest was the joint eigenvalue probability density function (JPDF) 
and the n-point correlation functions of GinUE, GinSE and GinOE. In the case of GinOE, 
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Ginibre was only able to report the JPDF in the restrictive case where all eigenvalues are 
real: One difficulty in determining the full JPDF being that the space of eigenvalues is 
naturally represented as a disjoint union indexed over the possible numbers of real and 
complex conjugate pairs of eigenvalues. The full JPDF for GinOE was finally computed in 
the 1990's by Lehmann and Sommers [6], and later independently by Edelman [2]. 

Surprisingly, the formulation for the ensemble average of a multiplicative class functions 
we present here is seemingly independent of this decomposition of space. Nonetheless, since 
we will use the JPDF (which is dependent on this decomposition) we introduce the details 
here. 

Throughout the discussion N will be fixed, and {L,M) will represent pairs of non- 
negative integers such that L + 2M = N. L will represent the number of real eigenvalues 
and M the number of complex conjugate pairs of eigenvalues. The symbols a and /3 will 
represent real and (non-real) complex eigenvalues respectively. We also use Xl and A2M to 
represent Lebesgue measure on and C*^ respectively. The space of all possible eigenvalues 
of N X N matrices can then be identified with the disjoint union 

y X (C \ R)^. 

{L,M) 

It shall be convenient to write the partial JPDF's in terms of the Vandermonde determinant. 
Given 7 e we define to be the N x N Vandermonde matrix in the coordinates of 7. 
(Superscripts on matrices will indicate variables on which the entries are dependent). The 
Vandermonde determinant is then given by 

(1.1) A(7):=dety^= n 7n-7m, 

Tn<n 

and given a € R^ and /3 G we define 

(1.2) A{a,/3) := detV^ where 7 := (/3i, ^1, . . . , /3m, /3m, ai, . . . , «/,) 
The partial JPDF is then given by Pl,m : R^ x (C \ R)^ ^ R where 

(1.3) P^,M(a,/3) = ej,' n n ^rf^ (v^|I^/3™l) e-('^^+^^)/^ 

e=i m=l 

and 

(1.4) e;v:=2^(^+i)/4nr(n/2). 

n=l 

This defines the full JPDF since the domains of the partial JPDF's are disjoint. The 
only essential difference between the formulation for P£, m(q:,/3) presented here and that 
presented by Lehmann and Sommers (and Edelman) is that they use the right hand side of 
(1.1). 

From (1.3) we can see the two main difficulties in the computation of (5'): (i) the 
decomposition of the space of eigenvalues, and (ii) the complicated nature of |A(a,/3)|. 
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2 Statement of Results 

In order to state our results it shall be convenient to define 

0(7) :=exp(-7V2) {erfc(V2|Im(7)|)}'^' ^(l)- 
We then define two skew-symmetric inner products, 

(P, Q)m. ■■= / <l){ai)<l){a2) P{ai)Q{a2) sgn(a2 - "i) dai da2, 

and 

(P,Q)c ■=-2i I PimiP) sgn(Im(/3))(iA2(/?) 



where sgn(O) — and sgn(a;) — x/\x\. These skew-symmetric inner products also appear in 
another paper on the statistics of the eigenvalues of matrices in GinOE by Kanzieper and 
Akemann [5]. 

Theorem 2.1. Let J be the integer part of {N + l)/2. //P = {Pi (7), Piil), Pn{i)} C 
C[7] is a set of monic polynomials with degP„ = n — 1, then, assuming (^') exists, 

(*) = e^i pfc/p, 

the Pfaffian of Up, where Up is the 2 J x 2 J antisymmetric matrix whose j,k entry is given 
by 

j {Pj,Pk)R + {Pj,Pk)c ifj,k<N, 
(2.1) Up[j,k] < sgn(A;-j) / (?i(a)Pmi„{j,fe}(a)rfa otherwise, 

and Cat is given as in (1-4)- 

Notice that when N is even the first condition in Equation (2.1) always holds. The 
asymmetry between even and odd cases is due to the fact that the Pfaffian is not defined 
for odd by odd matrices. 

If P satisfies the conditions of Theorem 2.1, we will say that P is a complete set of monic 
polynomials. When TV is even it is sensible to use a complete family of monic polynomials 
which are skew- orthogonal. 

Corollary 2.2. Suppose N = 2J, and let Q = {Qi, Q2, ■ ■ ■ ,Qn} be any complete family of 

monic polynomials specified by 

{Q2k-l,Q2j) = —{Q2j,Q2k~l) — Skj^j and {Q2j,Q2k) = {Q2j-1, Q2k-l) = 0, 

where (P, Q) = (P, Q)k + (P, Q)c- Then, 

J 

The quantities Tlj are referred to as the normalization(s) of Q. See [7, Ch. 5] or [3, 

Ch. 5] for more about skcw-orthogonal polynomials. 

If ip satisfies an additional symmetry, then we may write as a determinant 
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Corollary 2.3. Let J he the integer part of {N + l)/2. If ip{-l3) = for every /3 S C, 
and P is a complete family of monic polynomials in C[7] such that Pn is even when n — 1 
is even and odd when n — 1 is odd, then 

(*) = e^i dctAp 

where Ap is the J x J matrix whose j, k entry is given by 

Ap[j,k] :=C/p[2j-l,2fc]. 

We reiterate the striking fact that these formulations of (^1') are seemingly independent 
of the decomposition of the space of eigenvalues into all the different possible numbers of 
real and complex conjugate pairs of eigenvalues. In fact, the decomposition of the space of 
eigenvalues does enter into the statement of Theorem 2.1 — it is the reason that the inner 
product in the entries of Up are sums of skew inner products. One of these inner prod- 
ucts introduces to (vp) contributions from tp on real eigenvalues, while the other introduces 
contributions from pairs of complex conjugate eigenvalues. 

3 Averages over GinUE 

The simplicity of Theorem 2.1 and its corollaries can (perhaps) be most appreciated when 
compared with the analogous results for GinUE. GinUE is the set C^^^ together with the 
probability measure whose density is proportional to exp(Tr(Z*Z)/2). The joint eigenvalue 
probability density function was given by Ginibre as 

JV JV 

(3.1) PNh):=D],'\A{j)f]le-\'^"\"/\ where := {2nf l[rin + 1). 

n—l n—1 

We will denote the ensemble average of over GinUE by {^I/}, and we define the inner 
product, 

{P\Q) := / e-l^l'/2V(7)^Q(7)dA2(7). 

Theorem 3.1. Let P be any complete set of monic polynomials. Then, assuming {\E'} 
exists, 

W = D^i detl^p, 

where Wp is the N x N symmetric matrix whose j, k entry is given by Wp[j, k] = {Pj\Pk), 
and Djv is given as in (3.1). 

Corollary 3.2. Let Q be the complete family of monic polynomials which are orthogonal 
with respect to {■]■), then 

JV 

m = -D],^l[{Qn\Qn). 

n=l 

The proof of Theorem 3.1 is a straightforward exercise using well-known techniques of 
random matrix theory (apply the JPDF to the definition of {^}, expand A(7) as a sum 
over the symmetric groiip. apply Fubini's Theorem and simplify). A similar computation 
(in the domain of heights of polynomials) is done in [10]. The proof of 2.1 on the other hand 
is much more difficult due to the previous mentioned complications. The analogy between 
2.1 and 3.1 is made all the more striking when comparing the differences in difficulties of 
their proofs. Clearly there is a deeper phenomenon at work. 
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4 The Proof of Theorem 2.1 

The proof of Theorem 2.1 is based on, indeed almost identical to, another computation by the 
author in the study of heights of polynomials with integer coefficients [11]. The connection 
between that computation and the one presented here is that the Jacobian of the change of 
variables from the roots to the coefficients of a monic degree N real polynomial with L real 
roots and M pairs of complex conjugate roots is a constant times |A(a,/3)|/iy!A'/!. Since 
the audiences of these two results are likely disjoint many of the details are presented here, 
though the proof of many formulas which are purely combinatorial will only be referenced. 
From (1.3) we see, 

(4.1) (*) = 5^ / lf[4,{a)l[mm]\Ai<^,P)\dXL{a)dX2M{f3), 

Next, for each pair (L, M) we use the Laplace expansion in order to expand the Vandermonde 
determinant. Since the first 2M columns of V"'^ depend only on /3, while the remaining 
columns depend only on a, we will expand A(a, /3) via 2M x 2M and L x L minors. 



4.1 Notation for Minors 

For each K < N we define 3^ to be the set of increasing functions from {1,2,..., K} to 
{1,2..., TV}. That is, 

a^:={{l,2,...,K}^{l,2,...,7V} : t(l) < t(2) < ■ ■ ■ < t(K)}. 

Associated to each t G 3^ there exists a unique t' G 3]v-if such that the images of t and t' 
are disjoint. Each t € 3^ induces a unique permutation in G by specifying that 

i-iW I _ K <n<N. 

We define the sign of t by setting sgn(t) := sgn(it). The identity map in 3^ is denoted by i. 

Given an N x N matrix W and u, t G 3^, define Wu,t to be the K x K minor whose j, k 
entry is given by Wu,tb, ^1 = The complimentary minor is given by Wu',t'- 

As an example of the utility of this notation, the Laplace expansion of the determinant can 
be written as 

(4.2) det W = sgn(u) ^ sgn(l) det Wu,t • det Wu',t', 

where u is any fixed element of 3^. We will also use the abbreviated notation Wu for Wu,u', 
this is useful notation for working with Pfafiians since if W is an antisymmetric matrix then 
minors of the form Wu are also antisymmetric. 

We recall the definition of the Pfaffian here, li N = 2J and U is an N x N antisymmetric 
matrix, then the Pfaflaan of U is given by 



(4.3) 



Pff^=^ E sgn(r)n[/[T(2i-l),r(2i)], 
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where Sn is the Nth. symmetric group. We will also use the formula 

1 

(4.4) Pf^= J! E sgn(a)n Wi-l),T(2j)], 

■ o-enaj 3=1 

where II2J is the subset of S2J composed of those a with a{2j) > a{2j — 1) for j = 1, ... , J. 
The Pfaffian is related to the determinant by the formula det U = (Pf U)'^ (see for instance 
[9, Appendix: Pfaffians]). 

li U = R + C where R and C are antisymmetric 2 J x 2 J matrices, then Pf U has an 
expression in terms of the minors of R and C [11, Lemma 8.7]. 

J 

(4.5) Pf = E E • 

This will be useful since ultimately we intend to express (^) as the PfafBan of a sum of 
matrices. 

4.2 Steps in the Proof 

Lemma 4.1. Let 7 e C-'^ be given as in (1.2) and let W"'^ be the N x N matrix whose 
j, k entry is given by 

Then, if i G 3^ is the identity map on {1,2,..., 2M}, 
|A(a,/3)| = E sgn(t) LetW,^^{-i)^ JJ sgn9(/3„)| | det W^^t, [] sgn(afe - aA, 

where as suggested by the notation, the minors W^^ and W(^^, of W"'^ are dependent only 
on /3 and a respectively. 

Using Lemma 4.1, Equation (4.1) becomes 

W = e^' E ML! E / / [i{^{a,)\\\l<j^{f3M'^)\ 

{L,M) " "teaj'„ R^i^cM ^^=^ ^ '^'"=1 ^ 

x| det Wi?t, JJ sgn(afe - |detwft(-i)^ JJ sgnS>(/3„) I rfAL(a) rfA2M(/3), 

I j<fe J I m=l J 

and Fubini's Theorem yields 

(*) = E E / I n •^("^) n sg^^'^'^ - [ W^"t' dAi(a) 

(4.6) ^-^^ / |n de*^U^^2M(/3). 



CM 
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We remark that by using Fubini's Theorem we are assuming that the average {^) actually 
exists. 

Next, let K be the integer part of (L + 1)/2. Define T" to be the 2K x 2K antisymmetric 
matrix whose j, k entry is given by 

(4.7) T-[j,k]:={ i[ j,k<L + l, 
^ ' J 1^ sgn(K — otherwise. 

Then, it is well known that, 

(4.8) Yl sgn(afe-a,) =Pfr°, 

l<j<k<L 

(sec [3, Ch.5], or [11, Lemma 8.4]). It is worth remarking that when L is even then the first 
condition defining T" is always in force. Since the Pfaffian is only defined for even rank 
antisymmetric matrices, the second condition is used when L is odd to bootstrap a 2K x 2K 
antisymmetric matrix from an L x L matrix. 
Substituting (4.8) into (4.6) we see 

(*) = ^ J2 ^S^Wl! / |n'^("^)| Pf^" • detW^ft, dXL{a) 



Now, let J be the integer part of {N + l)/2. It is necessary for our calculations to replace 
the t G with elements of 32M- Each t G 3^ induces a unique to G t>y setting 
t = to. Notice that t' and differ in the fact that if TV ^ 2 J then t'o(2J - 2M) = 2J. 
Clearly, sgn(to) = sgn(t). 

Lemma 4.2. Let R be the 2 J x 2 J matrix whose j, k entry is given by 

{ {Pj,PkW if j,k<N + l 

-^b' ^] • \ ggn(/j _ f 0(q,) Pj^jjjjj (a) da otherwise, 

and suppose i G 3^ . Then, 

When A'^ is odd and t G then i?t/ is an odd by odd matrix. The introduction of to 
is useful since the Pfaffian of R^l is defined. 



Lemma 4.3. Let C be the 2J x 2 J matrix whose j, k entry is given by 
C\j,k] :-- 



{Pi,Pk)c if j,k<N + l 
otherwise. 
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and suppose i G 3^ . Then, 

' M 



M 



Ml 

c 



n <^(/3„)<^(/3„)sgn \ det dX2M{0) = Pf Ct„. 



L m=l 



Using Lemma 4.2 and Lemma 4.3 we may rewrite (4.9) as 

(4.10) (*) = ^ sgn(to)Pfi?t^ • PfCt^. 

{L,M) teJ^^ 

If u e then either 2 J is in the image of u or 2 J is in the image of u'. Notice that if 2 J 
is in the image of u then Pf Cu = 0. If 2 J is in the image of u' then u'(2 J — 2M) = 2J and 
hence u = to for some t e O^^. Thus we may replace the sum over 3^ in (4.10) with a 
sum over d'^M- Consequently, 

(*) = e^^ ^ J2 • 

(L,M)u63i^ 
J 

(4.11) = 51 E Pf • Pf 

where the second equation follows since the summand has been made to be independent of 
L. From (4.5) we see that (*) = e^^Pf(i? + C). Consequently, by (2.1), (*) = C^^ Pf C/p. 

4.3 The Proof of Lemma 4.1 

Applying (1.1) to (1.2), 



L M 



A{cx,(3) ^{Hiak- a,) > H 11 1^- " "'I 

j<k I ;=lm=l 



M 



And hence. 



(4.12) |A(a,/3)| = (-i)^ I nsgn(afc - a,) [] sgn3(/3„) i A(a,/3). 

I j<fe m=l J 

We may replace the monomials in the Vandcrmonde matrix with any complete family of 
monic polynomials without changing its dctcrmin ant. That is, A(a,/3) = dot VK""^. Using 
the Laplace expansion of the determinant (4.2) with u = i G , we see that 

detW«'^= J2 sgn(t)detWj°;'^-detWj?f. 
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Notice that the minors of the form W°['^ consists of elements from the first 2M columns of 
W"''^. These columns arc not dependent on a and thus wc may write these minors as W/'j. 
Similarly we may write the minors of the form W^ff as W^^j,. It follows that 

(4.13) detV°''^= sgn(t)detWi^j-detWi?t,, 

and the lemma follows by substituting (4.13) into (4.12) and simplifying. 

4.4 The Proof of Lemma 4.2 

We start with 

(4.14) 0= -i j^^detWi^^t' •Pfr°|n</)(a,)| d\L{cx), 

where t is an element of Expanding detW"i, as a sum over Sl allows us to write O 

as 

(4.15) O = -i 5] sgn(a) / \f[ct>{ae)f[Pt^e)ia^^ApfT"dXL{cx). 

^ V ' 

@ 

Sl naturally acts on by permuting the coordinates (denote this action by a • a), and it 
is easy to verify that for a € Sl, Pf T"^ " = sgn((T) Pf T". Using this fact, and the change 
of variables a ■ a we may write @ as 

@ = sgn(c7-i) / |n<^Mn^tw("^)l PiT'-dXLia). 
Substituting this into (4.15) we see that the sum over Sl exactly cancels 1/L\. That is, 

(4.16) 0= / TTL(a,)r[Ptm(a,)l PfT°dAL(a). 

Setting K to the integer part of {L + l)/2, and using (4.4), we see 

1 ^ 

(4.17) ¥iT'-=— J2 sgn(T)nsgn(a.(2fc)-a,(2fc_i)). 



Kl 

Ten2K k=i 



4.4.1 L Odd 



In the case of L odd, we view Ul+i = +oo so as to be consistent with (4.7). Substituting 

(4.17) into (4.16) we find 

1 f ^ ^ 

(4.18) ® = X] ^Sn(T) / JJ (j){ae)Pt(^e){ae) JJ sgn(Q!r(2fe) - a^(2k-i)) dXLia) . 
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For each r e 1X2^ there is a ko such that Q!^(2feo) = chl+i- If we set £0 = r(2fco — 1) then we 
may write © as 



K 



= </'(a^o)^t'(<o)("<?o)S n '/'("^)^t' Yl sgn(Q;^(2fe) - aT-(2fe-i)) 
t t=i k=i 



k=l 
k^ko 



X -P(t'or)(2/s)(ar(2/s))-P(t'or)(2fc-l)(Q!x(2fe-l)) Sgn(Q!^(2fe) - Q!r(2fe-1)) |, 

where the second equation comes from reindexing the first product by ^ i-^- T~^{t) together 
with the fact that 2{K — 1) = L — 1. Substituting this into (4.18) and applying Fubini's 
Theorem we find 

0=;^ sgn(T)^(^(x)P(t'0T)(2fe„-i)(a:)dx 

n / ^{^)^{y)P(i'oT)(2k){y)P(i'oT){2k-i){x)sg-B.{y - x)dxdy\ 

^ Sgn(T)|/" (l>{x)P(^i,oT){2k,-l){x)dx\ (-P(t'oT)(2fe-l),f(t'or)(2fc))M 



K 

X 



k^ko 

RecalUng the definition of gives (to o T)(2fco) = 2 J, and hence 

1 ^ 
(4.19) O = E SS^(^) ^t'o [^(2fco - 1), r(2fco)] [] [T(2fc - 1), T(2fc)] 

fc?=fco 

4.4.2 i Even 

When L is even, substituting (4.17) into (4.16) and simplifying, O becomes 



4t y] sgn(r)i TT / (/.(a;)^!>(?/)P(t/or)(2fc)(y)-P(i'or)(2fc-i)(a;)sgn(y-a;)da;dy 
^■rGn.^ \k=M J 

1 ^ 

(4-20) = Sgn(T) J|(P(t/or)(2fc-l),^'(t'or)(2fe))lR• 

' T-enas: fe=l 

Regardless if L is even or odd, (4.20) and (4.19) imply that, 

1 ^ 
®=:^ E sgn(r)n^t'„K2A:-l),r(2A;)]=Pfi?t- 
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4.5 The Proof of Lemma 4.3 

To prove Lemma 4.3 we set 



® = ^mT / i n 0(/3m)0(/3m)sgn3(/3„) I det dA2M(/3). 



From the definition of WP, we can write 



i,t 

M 



P{t0T){2m-l){Pm)P{i0T){2m){Pm)- 

TeS2M m=l 



Substituting this into ® we see 



® = E sgn(r)H)^^/ m 0(/J„)</>(/3„)sgn9(/3„)l 

-P(toT)(2m-l)(/3m)-P(to-r)(2m)(/3m) > dA2M(/3). 



By Pubini's Theorem, 

M 



® 2^M! 



■ TeS2M i. m=l '''^ 

xP{tor){2m-l)(P)P{tor){2m){l3) Sgn 9(/3)dA2 

1 ^ 
= E Sgn(T) ]J(P(tor)(2m-l),^(tor)(2m))c, 

" -r£S2M m=l 

which is Pf Ct- But, by definition, t = to, and hence ® = Pf Ct„ as desired. 



4.6 The Proof of Corollary 2.3 

Corollary 2.3 follows from the fact that if [/ is a 2J x 2J antisymmetric matrix, and 
A;] = when j — k = mod 2, then Pf ?7 = det A where A is the J x J matrix given by 
A[j, k] = U[2j — 1, 2k]. This is Lemma 8.10 of [11], or can be proven by conjugating U by 
an appropriate permutation matrix B so that 



A 
-A 



BUB-^ = 

Then, using well-known results about the PfafHan, PfU = P{{BUB-^) = det A. 
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